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In this paper, a new filtering method is presented to remove the Rician noise from magnetic resonance 
images (MRI) acquired using single coil MRI acquisition system. This filter is based on nonlocal neutro- 
sophic set (NLNS) approach of Wiener filtering. A neutrosophic set (NS), a part of neutrosophy theory, 
studies the origin, nature, and scope of neutralities, as well as their interactions with different ideational 
spectra. Now, we apply the neutrosophic set into image domain and define some concepts and operators 
for image denoising. First, the nonlocal mean is applied to the noisy MRI. The resultant image is trans- 
formed into NS domain, described using three membership sets: true (T), indeterminacy (I) and false 
(F). The entropy of the neutrosophic set is defined and employed to measure the indeterminacy. The co- 
Wiener filtering operation is used on T and F to decrease the set indeterminacy and to remove the noise. 
The experiments have been conducted on simulated MR images from Brainweb database and clinical 
MR images. The results show that the NLNS Wiener filter produces better denoising results in terms of 
qualitative and quantitative measures compared with other denoising methods, such as classical Wiener 
filter, the anisotropic diffusion filter, the total variation minimization and the nonlocal means filter. The 
visual and the diagnostic quality of the denoised image are well preserved. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Magnetic resonance imaging (MRI) is a powerful diagnostic 
technique used in radiology to visualize detailed internal struc- 
tures of the human body. MRI is primarily used to demonstrate 
pathological or other physiological alterations of living tissues and 
is a commonly used form of medical imaging 1]. Despite signif- 
icant improvements in recent years, MR images often suffer from 
low signal to noise ratio (SNR) especially in cardiac and brain imag- 
ing [2 . This is problematic for further tasks such as segmentation 
of important features; classification of images for computer aided 
diagnostics, three dimensional image reconstruction and image 
registration. Therefore, it is important to improve the SNR of the 
images used during quantitative analysis. This is usually done by 
using a denoising method, which is an important preprocessing 


Abbreviations: ADF, anisotropic diffusion filter; MRI, magnetic resonance imag- 
ing; NLM, nonlocal means; NLNS, nonlocal neutrosophic set; NS, neutrosophic set; 
PSNR, peak signal to noise ratio; SNR, signal to noise ratio; SSIM, structural similarity 
index; TV, total variation. 
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step used to improve the image quality by reducing the noise com- 
ponent while preserving all the image features. 

Numerous approaches of denoising MR images have been pro- 
posed in the literatures. Henkelman [3] was the first to estimate 
the magnitude MR image from a noisy image. McVeigh et al. [4] 
proposed the spatial filtering and temporal filtering for reduc- 
ing Gaussian noise in MR images. These filters eliminate the 
high frequency noise at the expense of blurring fine details and 
sharp edges in the MR images. A large number of edge-preserving 
methods have been proposed to overcome the above-mentioned 
blurring effects. For example, anisotropic diffusion filtering of 
MRI denoising proposed by Gerig et al. [5] are able to remove 
noise using gradient information while respecting important image 
structures. Recently, Krissian and Aja-Fernandez [6] proposed a 
new anisotropic diffusion filter based on a linear minimum mean 
squared error estimation and partial difference equations for Rician 
noise removal that has achieved state-of-the-art results. Samsonov 
and Johnson [7] proposed the noise adaptive non-linear diffusion 
method for filtering MR images with spatially varying noise levels. 
By using the receiver coil sensitivity profiles, the priori information 
regarding the image noise level spatial distribution is obtained and 
it is utilized for the local adjustment of the anisotropic diffusion 
filter. 

A simple wavelet based noise reduction was proposed by 
Weaver et al. [8]. The main drawback of the method was small 
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structures similar in size to noise were also eliminated. Wavelet 
transform based denoising for Rician noise removal in MRI was pro- 
posed by Nowak [9]. Pizurica et al. [10 proposed an alternative for 
denoising MR images with Rician noise using low complexity joint 
detection and estimation method. In the squared magnitude of the 
MR image, introduce an analytical model for the probability of sig- 
nal presence, which is adapted to the global coefficients histogram 
and to a local indicator of spatial activity (locally averaged magni- 
tude of the wavelet coefficients). Subtract the constant bias from 
the scaling coefficients, and subsequently compute the square root 
of the denoised squared-magnitude image. 

Bao and Zang [ 1 1 proposed a MR image denoising method based 
on an adaptive multiscale products threshold which incorporates 
the merits of interscale dependencies into the thresholding tech- 
nique for denoising. Two adjacent wavelet subbands are multiplied 
to amplify the significant features and to dilute the noise. Apply 
thresholding to the multiscale products instead of the wavelet 
coefficients to better differentiate edge structures for noise. The 
distribution of the products analyzed and an adaptive threshold is 
formulated to remove most of the noise. Thus, it is a good edge pre- 
serving denoising method. Yang and Fei [12] proposed the wavelet 
multiscale denoising method for the Rician nature of MR data. In 
this method, Radon transform is applied to the original MR images 
and the Gaussian noise model is used to process the MR sinogram 
image. A translation invariant wavelet transform is employed to 
decompose the MR sinogram into multiscales in order to effectively 
denoise the image. 

Nonlocal means (NLM) filter introduced by Buades et al. [13] 
has emerged as a very simple and effective way to reduce noise 
while minimally affecting the original structures of the image. Man- 
jon et al. [14] modified the original NLM algorithm to denoise 
the multispectral MRI. In the multispectral sequences, the similar- 
ity measure can be obtained by combining information of various 
channels. Manjon et al. [15] proposed the unbiased NLM approach 
for MRI denoising. The unbiased NLM is obtained by subtracting 
the noise bias from the squared value of NLM. The main drawback 
of the NLM algorithm is computational burden due to its complex- 
ity especially on 3D MRI data. In order to overcome this, Coupe 
et al. [16] proposed an optimized blockwise NLM filter for denois- 
ing 3D MRI. This is achieved by the following steps: an automatic 
tuning of the smoothing parameter, a selection of the most relevant 
voxels for NL means computation, a blockwise implementation 
and a parallelized computation. Manjon et al. [17] proposed the 
denoising methods for three-dimensional MR images by exploiting 
the sparseness and self-similarity properties of the images. These 
methods are based on a three-dimensional moving window cosine 
transform hard thresholding and a three-dimensional rotationally 
invariant version of the NLM filter. 

Sijbers et al. [18-20] estimated the Rician noise level and 
performed signal reconstruction using maximum likelihood (ML) 
approach for reducing bias that appears in the conventional 
approach. Sijbers et al. [21 used this approach to estimate the 
image noise variance from the background mode of the histogram 
of MR image which is known to be Rayleigh distributed. He and 
Greenshields [22] used the nonlocal maximum likelihood (NLML) 
estimation method for Rician noise. This method is based on 
deploying maximum likelihood estimator on the nonlocal neigh- 
borhood in order to predict the underlying noise. Rajan et al. 
[23 used the ML based estimation of the local variance for each 
pixel of the image using a local neighborhood when no back- 
ground information is available, like cardiac and lung images. 
Rajan et al. [24] proposed the MRI denoising for spatially vary- 
ing noise levels based on ML estimation using restricted local 
neighborhoods. Recently, Rajan et al. [25] presented nonlocal 
maximum likelihood estimation for denoising multiple coil MR 
images. 


This paper focuses on developing a new denoising method based 
on the nonlocal neutrosophic set (NLNS) approach of Wiener filter- 
ing to remove the Rician noise in MRI. First the NLM method is 
applied to the noisy MRI and the resultant image is transformed 
into the neutrosophic set. The co-Wiener filtering is employed to 
reduce the indetermination degree of the image, which is measured 
by the entropy of the indeterminate subset, after filtering, the noise 
is removed. The rest of this paper is organized as follows: Section 
2 discusses the characteristics of the noise in MRI. In Section 3, 
the proposed methodology of MRI denoising is described in detail. 
The materials used to evaluate the performance of the proposed 
method, the validation strategies and the comparative evaluation of 
the denoising results are discussed in Section 4. Section 5 presents 
the conclusion and the future scope. 

2. Characteristics of the noise in MRI 

The raw data obtained during MRI scanning are complex values 
that represent the Fourier transform of a magnetization distribu- 
tion of a volume of tissue at a certain point of time. An inverse 
Fourier transform converts these raw data into magnitude, fre- 
quency and phase components that more directly represent the 
physiological and morphological features of interest in the person 
being scanned. So, noise in the k-space in MR data from each coil 
is assumed to be a zero mean uncorrelated Gaussian process with 
equal variance in both real and imaginary parts because of the lin- 
earity and orthogonality of the Fourier transform. However, it is 
common practice to transform the complex valued images into 
magnitude and phase images. Since computation of a magnitude 
(or phase) image is a nonlinear operation, the probability density 
function (PDF) of the data under concern changes. The magnitude 
data in spatial domain is modeled as the Rician distribution and 
the so-called Rician noise (the error between the underlying image 
intensities and the measurement data) is locally signal dependent 
[26]. 

The complex MR data is given as, 

M = a + ib (1) 

The noise added to the complex raw data is zero mean Gaussian 
noise N( 0, a 2 ), 

M = (a + n re ) + i(b + Him) (2) 


The spatial MR image is the magnitude of the noisy raw data. There- 
fore, 



Thus, the distribution of |M| becomes Rician [23] and is represented 
as, 


Pm(M \a, o n ) = ^e-W 2+A2 V 2a h 0 ( ^f) u(M ) ( 4 ) 

' °n \ a n J 

with Jo(-) as the modified zeroth order Bessel function of the first 
kind, o 2 is the noise variance, A is the noiseless signal level, M is 
the MR magnitude variable and u(-) the Heaviside step function. 
In high SNR, i.e., high intensity (bright) regions of the magnitude 
image, the Rician distribution tends to a Gaussian distribution with 
mean a/a 2 + o 2 and variance a 2 given as 

Pm(M \A, o n ) « 2 e ~ (iv|2 - \A^+^V 2g n u(M) ( 5 ) 

v 27ror n 
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In the image background, where SNR is zero due to the lack of water 
proton density in the air, the Rician PDF simplifies to a Rayleigh 
distribution with PDF 

Pm(M, o„) = ^ e -(M 2 /2<r„ 2 ) u(M) (6) 

°n 

Thus, Rician noise in MR images behaves to be Gaussian distributed 
when SNR is high and Rayleigh distributed for low SNR. The squared 
magnitude image (the value of each pixel in the image is the square 
of the value of the corresponding pixel in the original magnitude 
image) has a noise bias which is signal independent and it can be 
easily removed. Such bias is equal to 2a 2 as shown by Nowak [9] 
and therefore, a simple bias subtraction will recover its original 
value. The estimation of the noise standard deviation o from the 
noisy MR image is based on the local estimation of the skewness of 
the magnitude data distribution proposed by Rajan et al. [23]. 


II 1 1 2 

weighted Euclidean distance, ||u(Nj) - u(Nj) || , where a>0 is the 

standard deviation of the Gaussian kernel. The application of the 
Euclidean distance to the noisy neighborhoods raises the following 
equality 

E||u(N0-u(N,.)|£ a = ||u(N i )-u(N J )||^ + 2a 2 (8) 

This equality shows the robustness of the algorithm since in expec- 
tation the Euclidean distance conserves the order of similarity 
between pixels. 

The pixels with a similar gray level neighborhood to u(N z ) have 
larger weights in the average. These weights are defined as, 





|u(N i )-u(N J -)|| / h 2 


(9) 


where Z(i) is the normalizing constant. 


3. Proposed methods 

Neutrosophy, a branch of philosophy introduced in [27] as a 
generalization of dialectics, studies the origin, nature and scope of 
neutralities, as well as their interactions with different ideational 
spectra. In neutrosophy theory, every event has not only a certain 
degree of the truth, but also a falsity degree and an indeterminacy 
degree that have to be considered independently from each other 
[27]. Thus, a theory, event, concept, or entity, {5} is considered with 
its opposite {Anti-5} and the neutrality {Neut-S}. {Neut-5} is nei- 
ther {5} nor {Anti-5}. The {Neut-5} and {Anti-5} are referred to 
as {Non-5}. According to this theory, every idea {5} is neutralized 
and balanced by {Anti-5} and {Non-5} [27]. Neutrosophic set (NS) 
provides a powerful tool to deal with indeterminacy. 

NS had been applied to image thresholding, image segmenta- 
tion and image denoising applications [28-32]. Cheng and Guo [28] 
proposed a thresholding algorithm based on neutrosophy, which 
could select the thresholds automatically and effectively. The NS 
approach of image segmentation for real images [29 and for color 
texture image segmentation discussed in [30]. In [31], some con- 
cepts and operators were defined based on NS and applied for image 
denoising. It can process not only noisy images with different lev- 
els of noise, but also images with different kinds of noise well. 
The method proposed in [31 ] adapted and applied for MRI denois- 
ing [32], which shows some significant improvement in denoising. 
This paper is presented as an extended work of [33 ] ; MRI denoising 
based on NLNS based Wiener filtering is evaluated with Simulated 
and clinical MR images. Since the NLM filtered image is used as the 
reference image for NS domain to reduce the indeterminacy of the 
sets, this method is suitable for denoising MR images with Rician 
noise. 

3.1. Nonlocal means operation 

The NLM method is applied to the noisy MRI in order to create 
the reference image. In the nonlocal means [13], Given a discrete 
noisy image u = (u(i) | i e /}, the estimated value NL [u] (i), for a pixel 
i, is computed as a weighted average of all the pixels in the image, 


rj, .x -||“(Ni)-u(Nj)||^ /h 2 

Z(i) = > e 11 ll2 >« (10) 

j 

and the parameter h acts as a degree of filtering. It controls the 
decay of the exponential function and therefore the decay of the 
weights as a function of the Euclidean distances. 

3.2. The image in neutrosophic set 

The NLM filtered image is transformed into NS domain. A neu- 
trosophic set image P NS is characterized by three membership sets 
T,/,F. A pixel P in the image is described as P(T,/,F) and belongs to 
W in the following way: it is t true in the set, i indeterminate in 
the set, and /false in the set, where t varies in T, i varies in I and / 
varies in F. Then the pixel P(ij) in the image domain is transformed 
into the neutrosophic set domain PNs(hj) = {T{i,j), /(z, j), F(z, j)}. T(i, 
j), J(i, j) and F(z, j) are the probabilities belong to white pixels set, 
indeterminate set and non white pixels set, respectively [31-33], 
which are defined as: 


T(i, j) = M f min (11) 

gmax-gmin 

i+w/2 j+w/2 

g(U) = — 1— V' V' g(m, n) (12) 

w x W ^ — s z — s 

m=i-w/2n=j-w/2 

HU) = f J) ~f min (13) 

d max — min 

S(i,j) = abs(g(i,j)-g(i,j )) (14) 

F(i,j)=\-T(i,j ) (15) 


where g(i,j) is the local mean value of the pixels of the window. S(i, 
j) is the absolute value of difference between intensity g(z, j) and its 
local mean value g(z, j). 


NL[u](i) = y W (iJ)u(J) (7) 

hi 

where the family of weights {w(z,j)}j depends on the similar- 
ity between the pixels i and j, and satisfy the usual conditions 
0 < w(i, j) < 1 and ^\w(i,j) = 1. 

The similarity between two pixels i and j depends on the sim- 
ilarity of the intensity gray level vectors u(N 2 ) and iz(Nj), where 
N k denotes a square neighborhood of fixed size and centered at 
a pixel k. This similarity is measured as a decreasing function of the 


3.3. Neutrosophic image entropy 

For an image, the entropy is utilized to evaluate the distribution 
of the gray levels. If the entropy is maxim, the intensities have equal 
probability. If the entropy is small, the intensity distribution is non- 
uniform. 

Neutrosophic entropy of an image is defined as the summation 
of the entropies of three subsets T, I and F [3 1 -33 ] : 

Fn NS = Enj + Frz/ + Frz f 


(16) 
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max{F} 


En T = - y ] p T (a)lnp T (a) 

(17) 

a=min{T} 


max{/} 


En, = - y ] p,(a)lnp,(a) 

(18) 

a=min{/} 


max{F} 


En F = - y ] Pf(a)lnp F (a) 

(19) 

a=min{F} 


where En T , Exi\ and En F are the entropies of sets 

T, / and F, respec- 

tively. p F (a ), p/(a) and p F {a) are the probabilities 
and F, respectively. 

of elements in T, I 

3.4. c o-Wiener filtering operation 


The values of /(z, j) are employed to measure the indeterminate 
degree of element Pjvs(h j)- To make the set I correlated with T and 
F, the changes in T and F influence the distribution of element in 
I and vary the entropy of I. A co-Wiener filtering operation for P NS , 

P NS (col is defined as [33]: 
P N s(co) = P(t(co),I(colP(co)) 

(20) 

{ T I < co 

n<o) = 

[ T&) I > co 

(21) 

Tco(iJ) = wiener [T(m, n)} 

(m,n)eSjj 

(22) 

[F I < CO 

* r 

[^ Fqj I > co 

(23) 

F M (i,j) = wiener {F(m, rz)} 

(m,n)eS ij 

(24) 

UiJ)= S f J) ~ S s tmin 

T max T min 

(25) 

8 t (i,j) = abs(t(iJ)-t(i,j )) 

i+w/2 j+w/2 

t(U)=r-2 L- V' V f(m, n) 
w x W 

(26) 

(27) 

m=i-w/2n=j-w/2 



where is the absolute value of difference between intensity 

f (z, j) and its local mean value T(z, j) at (z, j) after co-Wiener filtering 
operation. 

After new concepts and operators are defined in the neutro- 
sophic set for an image, a MRI denoising method based on nonlocal 
neutrosophic set approach of Wiener filtering is proposed and sum- 
marized as below (see Fig. 1): 

(1 ) NLM algorithm is applied to the noisy MRI; 

(2) Transform the resultant image into NS domain; 

(3) Use co-Wiener filtering operation on the true subset T to obtain 

T 'col 

(4) Compute the entropy of the indeterminate subset J w , En? (z); 

1(0 

(5) if {{En, (i + 1 ) - En-. ( i))/En-. (i)) < S, go to Step 5; else T = f a , 

1(0 1(0 1(0 

go to Step 2; 

(6) Transform subset T M from the neutrosophic domain into the 
gray level domain. 


Noisy image 



Apply Non Local Mean Algorithm. 


Transform the NLM filtered image 
into Neutrosophic Set Domain 


Apply the wiener filter to reduce 
the indeterminacy 





Transform NS domain to gray level 

domain 


Fig. 1 . Block diagram of the MRI denoising method based on NLNS Wiener filtering. 

4. Experimental results and discussions 

4.1. Materials 

The experiments were conducted on two MRI datasets. The first 
data set consists of simulated MR images obtained from the Brain- 
web database [34]. The second data set consists of clinical MRI 
collected from PSG Institute of Medical Sciences and Research (PSG 
IMS & R), Coimbatore, Tamilnadu, India. Simulated MR images are 
used collectively as the reference to evaluate and compare the 
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fry 


ADF Residual 



Wiener Residual 



TV Residual 



Fig. 2. Denoising results and the corresponding absolute value of residuals for simulated T1 weighted axial MRI corrupted by 7% of Rician noise. 


validity of the proposed technique. It evicts the data dependency 
enabling precise comparative studies. The data set consists of T1 
weighted Axial, T2 weighted Axial and T1 weighted Axial with 
multiple sclerosis (MS) lesion volumes of 181 x 217 x 181 voxels 
(voxels resolution is 1 mm 3 ), which are corrupted with different 
level of Rician noise (1% (a = 2.5) to 15% (a = 37.5) of maximum 


intensity). Rician noise was generated by Gaussian noise to real and 
imaginary parts and then computing the magnitude of the image. 

In the clinical data sets, the images were acquired using Siemens 
Magnetom Avanto 1 .5 T Scanner and Hitachi 0.3 T AIRIS II scanner. 
MRI pulse sequences are pattern of radiofrequency pulses and mag- 
netic gradients that are used to produce an image. Here, the images 
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Original Noisy 



ADF Wiener TV 





ADF Residual 


Wiener Residual 



NLM NLNS 


TV Residual 





NLM Residual 


NLNS Residual 




Fig. 3. Denoising results and the corresponding absolute value of residuals for simulated T2 weighted axial MRI corrupted by 9% of Rician noise. 


are acquired using spin echo sequences with long repetition time 
(TR) and short echo time (TE). 

1. Siemens Magnetom Avanto 1.5 T Scanner 

a. T2 weighted Sagittal MR image of normal brain with 
TR = 4460ms,TE = 85ms, 5 mm thickness and 512 x 512 reso- 
lution. 


b. T1 weighted Axial MR image of tumor pathology with 
TR = 550ms, TE = 8.7ms, 5 mm thickness and 499x559 res- 
olution. 

2. Hitachi 0.3 T AIRIS II Scanner 

a. T1 weighted Axial MR image of tumor pathology with 
TR = 249 ms, TE = 1 8 ms, 6.5 mm thickness and 256x256 res- 
olution. 
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TV Residual 








NLM Residual NLNS Residual 




Fig. 4. Denoising results and the corresponding absolute value of residuals for simulated T1 weighted axial MRI with MS lesion corrupted by 15% of Rician noise. 


b. T2 weighted Coronal MR image of normal brain with 
TR = 5472 ms, TE = 1 25 ms, 7 mm thickness and 256x256x 18 
resolution. 

4.2. Validation strategies 


the structural similarity (SSIM) index [35]. The peak signal to noise 
ratio in decibel (dB) is measured using the following formula: 


PSNR = -10 log 10 


E i=H — 1 

i = 0 


H xWx 255 2 


(28) 


The performance of the denoising algorithm is measured by 
using the quality metrics such as peak-signal-to-noise ratio (PSNR), 


where J(i, j) and /^(i, j) represent the intensities of pixels (i, j) in 
the original image and denoised image, respectively. The higher 
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Tlw 



% of Noise level 



T2w T2w 




TlwMS TlwMS 




Fig. 5. Comparison results for Brainweb simulated MR images. (Left) PSNR of the compared methods for different image types and noise levels. (Right) SSIM of the compared 
methods for different image types and noise levels. 


the PSNR, the better the denoising algorithm is. SSIM gives the 
measure of the structural similarity between the original and the 
denoised images and are in the range of 0 to 1. The SSIM works 
as follows: let x and y be two non negative images, where as one 
has perfect quality. Then, the SSIM can serve as a quantitative mea- 
sure of the similarity of the second image. The system separates the 
task of similarity measurement into three comparisons: luminance, 
contrast and structure. It can be defined as 


SSIM(x,y) 


(2/Xx/Xy + Cl )(2<7 xy + C 2 ) 
(Mx + My + Cl )(cr^ + CTy + C 2 ) 


(29) 


where /jl x and /x y are the estimated mean intensity and o x and o y 
are the standard deviations, respectively. o xy can be estimated as 

1 N 

a xy = 1 — MxXyi — My) (30) 

i=l 

Ci and C 2 in Eq. (29) are constants and the values are given as 
Ci = (/Ci I) 2 and C 2 = (/C 2 L) 2 where K \ , I< 2 < 1 is a small constant and 
L is the dynamic range of the pixel values (255 for 8 bit gray scale 
images). 
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Original noisy ADF Wiener 



Fig. 6. Denoising results of clinical T2 weighted Sagittal MRI. 


The residual image is also obtained by subtracting the denoised 
image from the noisy image [15]. The residual image is required to 
verify the traces of anatomical information removed during denois- 
ing. Hence, this reveals the excessive smoothing and the blurring 
of small structural details contained in the image. 

4.3. Performance comparison 

The performance of the proposed method has been com- 
pared with several standard and established algorithms such as 


Wiener filter, anisotropic diffusion filter, total variation minimi- 
zation scheme and nonlocal means filter. 


4.3 A. Wiener filter 

N. Wiener proposed the concept of Wiener filtering also 
called as minimum mean square error filter or least square 
error filter in 1942 36]. This method is founded on consider- 
ing images and noise as random processes and the objective is 
to estimate / of the uncorrupted image / such that the mean 



Fig. 7. Denoising results of small part of the clinical T1 weighted axial MRI with tumor pathology. 












788 


J. Mohan et al. / Biomedical Signal Processing and Control 8 (2013) 779-791 


T1 weighted axial Brainweb database Volume of images 




23 



16 1 1 1 1 1 

40 50 60 70 80 90 100 

MRI Slice Number 
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Fig. 8. Results of experiments performed on a sequence of MR images collected from Brainweb MR volume. 


square error between them is minimized. This error measure is 
given by 

e 2 =E{(f-f) 2 } (31) 


where F {•} is the expected value of the argument. In the frequency 
domain, 


F(u, v) = 


H*(u, v)Sf{u, v ) 

Sf{u, i/)| H(u, v)\ 2 +S v {u, v) 


G(u, v) 


(32) 


where H(u,v) is the transform degraded function, H*(u,y) is the com- 
plex conjugate of H(u,v\ S n (u, v) is the power spectrum of the noise, 
Sf(u,v) is the power spectrum of the undegraded image and G(u,v) 
is the transform of the degraded image. The restored image in the 


spatial domain is given by the inverse transform of the frequency 
domain estimate F(u, v). 


4.3.2. Anisotropic diffusion filter (ADF) 

Perona and Malik [37] introduced the anisotropic diffusion fil- 
ter. In this approach the image u is only convolved in the direction 
orthogonal to the gradient of the image which ensures the preser- 
vation of edges. The iterative denoising process of initial image u 0 
can be expressed as 

f 9U Q X f ’ - = div(c(x, t)Vu(x, t)) 

I u(x, 0) = u 0 (x) 
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Table 1 

Comparison of the denoising techniques based on the performance metrics for simulated MR images. 
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MR image 

Denoising methods 

Performance metrics 
PSNR (dB) 

SSIM 

T1 -weighted brain MRI corrupted by 7% Rician noise 

Wiener 

24.4 

0.9429 


TV 

24.5 

0.9540 


ADF 

25.3 

0.9643 


NLM 

25.32 

0.9646 


NLNS 

31.35 

0.9894 

T2-weighted brain MRI corrupted by 9% Rician noise 

Wiener 

22.02 

0.8334 


TV 

21.58 

0.9342 


ADF 

22.64 

0.9468 


NLM 

23.37 

0.9494 


NLNS 

27.11 

0.9813 

T1 -weighted brain MRI with MS lesion corrupted by 

Wiener 

17.07 

0.7870 

15% Rician noise 

TV 

16.9 

0.7901 


ADF 

16.84 

0.8127 


NLM 

16.88 

0.8370 


NLNS 

23.92 

0.9254 


where Vu(x, t) is the image gradient at voxel x and iteration t, du(x, parameters were chosen as w = 5, /= 2 and h is proportional to the 

t)/3t is the partial temporal derivation of u(x, t) and noise level of the image [13]. 


c(x, t) = g(|| Vu(x, t)||) = e II Vu (*4)||/k 2 ( 34) 

where I< is the diffusivity parameter. 


4.3.3. Total variation (TV) minimization scheme 

The difficult task to preserve edges while correctly denoising 
constant areas has been addressed also by Rudin, Osher and Fatemi. 
They proposed to minimize the TV norm subject to noise con- 
straints [38], that is 

u = argmin / |Vu(x, t)\ dx (35) 

uei2 3 J 


subject to 

/ (u(x) - u 0 (x))dx = 0 and / |u(x) - u 0 (x)| 2 dx = o 2 (36) 

Jn 3 Jqs 

where u 0 is the original noisy image, u is the restored image and a 
is the standard deviation of the noise. In this model, the TV minimi- 
zation tends to smooth inside the image structures while keeping 
the integrity of boundaries. The TV minimization scheme can be 
expressed as an unconstrained problem 


u = argmin 

uei2 3 


J | Vu(x)| dx + A J |u(x)-u 0 (x)| 2 dx (37) 


where A is a Lagrange multiplier which controls the balance 
between the TV norm and the fidelity term. Thus A acts as the 
filtering parameter. Indeed, for high values for A the fidelity 
term is encouraged. For small values for A the regularity term is 
desired. 

All the methods were implemented using MATLAB 2010a (The 
Math Works, Inc.) in Windows XP 64-bit Edition (Pentium Dual Core 
2.4GFIz with 4 GB of RAM). Wiener filtering was computed using 
the function from Matlab Image Processing Toolbox [39]. For AD fil- 
ter, the parameter I< varies from 0.05 to 1 with a step of 0.05 and the 
number of iterations varies from 1 to 15. For TV minimization, the 
parameter A varies from 0.01 to 1 with step of 0.01 and the num- 
ber of iterations varies from 1 to 1 0. There are three key parameters 
that need to be set to use the NLM filter; one is the size of the search 
window (w), second one is the size of the neighborhood window 
(f) and finally, the degree of the filtering (h). For the simulation, the 


4.4. Evaluation on simulated data set 

The detailed images of the denoising results obtained for the 
T1 weighted, T2 weighted and T1 weighted with MS lesion axial 
images corrupted by 7% (a = 17.36, PSNR = 20.52 dB), 9% (a = 22.32, 
PSNR= 18.31 dB) and 15% (or = 37.5, PSNR= 13.88 dB) of Rician noise 
respectively shown in Figs. 2-4. From the visual comparison of 
these results, the proposed method surpassed all other methods 
at high noise levels on the three types of data in terms of producing 
more detailed denoised image in which all the distinct features and 
small structural details are well preserved. 

The PSNR and the SSIM values obtained for the T1 weighted, 
T2 weighted and T1 weighted with MS lesion axial images with 
different noise levels ( 1 % (or = 2.5, PSNR = 37.39 dB) to 1 5% (or = 37.5, 
PSNR= 13.88 dB)) using the aforementioned denoising techniques 
are given in Fig. 5. As the level of noise increase, the performance of 
the proposed NLNS Wiener filter shows significant improvement 
over the other denoising methods. Table 1 shows a comparison of 
the experimental results for the denoising methods based on PSNR 
and SSIM for the T1 weighted, T2 weighted and T1 weighted with 
MS lesion axial images corrupted by 7%, 9% and 15%, respectively. 
Fligher the value of PSNR and higher the value of SSIM shows 
that the proposed filter perform superior than the other denoising 
methods. 

4.5. Evaluation on clinical data set 

The detailed images of the denoising results obtained for the 
T2 weighted Sagittal MR image of normal brain are given in Fig. 6. 
The detailed features and edges in the image are well preserved 
by the proposed method compared with other denoising methods. 
The denoising results obtained for the small part of T1 weighted 
axial brain with tumor pathology are shown in Fig. 7. With the 
other denoising methods (ADF, Wiener, TV and NLM), the sharp- 
ness along the edges and the pathology structure is extensively 
blurred. Comparing the NLM and NLNS Wiener, it is evident that 
NLNS Wiener preserves well the details. 

4.6. Evaluation on sequence of images 

The performance was also evaluated for the sequence volume of 
MRI scan at the specific intervals in order to take into account the 
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influence of structural variations of the efficiency of the denoising 
techniques. The results for the simulated MRI volumes of the T1 
weighted slices of 40-100 with 7%, 9% and 1 5% noise levels, respec- 
tively, shown in Fig. 8. It can be noted that the performance of the 
proposed method is superior to the other denoising techniques. 


5. Conclusion 

In this paper, a novel MRI denoising technique based on nonlocal 
neutrosophic set approach of Wiener filtering has been proposed. 
The method is adapted to Rician noise characteristics. The NLM 
filtered MR image is described as a NS set using three membership 
sets T, J and F. The entropy in neutrosophic image domain is defined 
and employed to measure the indetermination. The Wiener filter 
is applied to reduce the set’s indetermination and to remove the 
noise in the MR image. 

The performance of the proposed denoising filter is compared 
with Wiener, ADF, TV and NLM based on PSNR and SSIM. The 
experimental results demonstrate that the proposed approach can 
remove noise automatically and effectively. This filtering method 
tends to produce good denoised image not only in terms of visual 
perception but also in terms of the quality metrics. In the clinical 
MRI with tumor pathology, the filter preserves the major visual 
signature of the given pathology. Even though the computational 
complexity of the proposed method is equal to the NLM filter, on the 
contrary it increases PSNR and SSIM value of the denoised image. 
And finally, further works should be implementing spatially adap- 
tive Wiener filtering in the neutrosophic domain for accounting the 
spatially varying noise level in parallel imaging MRI. 
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